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We show that pristine M 0 S 2 single layer (SL) exhibits two bandgaps Eg\\ =1.9 eV and Egi_ = 3.2 
eV for the optical in-plane and out-of-plane susceptibilities X|| and respectively. In particular, we 
show that odd states bound to vacancy defects (VDs) lead to resonances in x± inside Eg± in M 0 S 2 
SL with VDs. We use density functional theory, the tight-binding model, and the Dirac equation to 
study M 0 S 2 SL with three types of VDs: (i) Mo-vacancy, (ii) S 2 -vacancy, and (iii) 3xMoS2 quantum 
antidot. The resulting optical spectra identify and characterize the VDs. 


PACS numbers: 61.72.jd,42.65.An,71.15.-m,73.22.-f 

Introduction. Monolayer transition metal dichalco- 
genides (TMDCs) (MX2; M= transition metal such as 
Mo, W and X=S, Se, Te) have attracted a lot of atten¬ 
tion due to their intriguing electronic properties. Mono- 
layer TMDCs are semiconductors with direct bandgap 
Eg\\ in the visible range, which makes them suitable for 
optoelectronic, spintronic, valleytronic, and photodetec¬ 
tor devices In order to increase the performance 

of such devices based on TMDC single layer (SL), it is 
crucial to characterize the defects present in TMDC SLs. 

Here we show that the bandgap Eg± = 3.2 eV for the 
optical out-of-plane susceptibility xa pristine M0S2 
SL provides a large energy window to characterize va¬ 
cancy defects (VDs). Pristine M0S2 SL is invariant with 
respect to ay reflection about the z = 0 (Mo) plane, 
where the a axis is oriented perpendicular to the Mo 
plane. Therefore, electron states break down into two 
classes: even and odd, or symmetric and anti-symmetric 
with respect to ay. We show below that this leads to 
the nontrivial consequence that X;s = Xa has a bandgap 
of Egj_ = 3.2 eV, which is substantially larger than the 
bandgap Eg\\ = 1.9 eV for the in-plane component of 
the optical susceptibility Xx = Xy = X||- we show, 
due to the optical selection rules for the even and odd 
states, there are no tt transitions, driven by A-polarized 
photons, below 3.2 eV. Hence, xa foi* pristine M0S2 SL 
must vanish for energies below 3.2 eV. 

Several studies on VDs in 2D materials have emerged. 
The minibands resulting from quantum antidot (QAD) 
superlattices can be used to tune the bandgaps of 
graphene [8] and M0S2 SL Ellin]. In another study, 
we have shown that substitutional defects in the form of 
M0O3 not only lead to strong suppression of the conduc¬ 
tivity El but also to photoluminescence quenching [12]. 
Recently, VDs in M0S2 SL have been characterized the¬ 
oretically in terms of magnetic properties m- A recent 
experimental study used scanning transmission electron 
spectroscopy to characterize several types of defects in 
M0S2 SL, including Mo, S, and S2 VDs [Il|. M0S2 SL 
with S-vacancies might catalyze alcohol synthesis from 


syngas US). 

Here we show that VDs yield strong resonances in xa , 
which provides the opportunity to optically characterize 
VDs in M0S2 SL with VDs (denoted by M0S2 SLVD). 
We consider the optical signatures of states bound to 
three types of VDs in M0S2 SLVD: (i) Mo-vacancy, (ii) 
S2-vacancy, and (iii) a hexagonal 3xMoS2 QAD (see Fig. 

0 . 

Bandstructure. First we start with the numerical 
bandstructure calculation of M0S2 SLVD using standard 
Density Functional Theory (DFT) with meta-GGA func¬ 
tionals m, providing accurate estimates of bandgaps 
without the need to perform computationally inten¬ 
sive DFT calculations using the GW approximation 
iniiiH]. The calculations are implemented within Atom- 
istix Toolkit 2014.2 [ 19 ]. The resulting bandstructures 
are shown in Fig. The periodic structure of the su¬ 
perlattice allows one to characterize the electron states 
by the bandstructure (k) , where k is the vector in the 
first Brillouin zone of the superlattice and n enumerates 
different bands. We consider supercells with dimensions 
7 x 7 x 1 (Fig.| b, c) and 8x8x1 (Fig. d) having 
147 and 192 number of atoms, respectively. For Brillouin 
zone integration we consider k sampling of 7 x 7 x 1. 
The cut off energy is set to 300 eV and the structure is 
optimized by using a force convergence of 0.01 eV/A. 

Tight-binding model (TBM) and symmetries. Within 
the TBM approximation the electron wavefunction can 
be presented as where 

j enumerates atoms in the layer and the summation over 
/i runs over respective atomic orbitals, whose set for the 
j-th atom is denoted Oj. Choosing Bx^y in the plane 
of the layer and Bz perpendicularly, for Mo the real or¬ 
bitals of main importance are the d-orbitals (7^,2 _^2, dxy 
and so on, while for S atoms these are p-orbitals 
with i = x,y,z and t and b denoting the top and bot¬ 
tom layers, respectively. The classification of the elec¬ 
tron states simplifies when the symmetry with respect 
to ay : z —z is taken into account. The electron 
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FIG. 1. (a) Bandstructure of M 0 S 2 SL, showing the out-of-plane bandgap Eg± — 3.2 eV determined by the transition Tq. The 
Fermi level is set at cf = 0 eV. (b) Mo-vacancy in 7x7 supercell, (c) S 2 -vacancy consisting of a pair of S atoms removed in 7x7 
supercell, (d) hexagonal 3xMoS2 QAD in 8x8 supercell. 


states transform according to A and M', the irreducible 
representations of Z2 = {e^ah}- The respective even 
and odd orbitals are locally spanned by the bases [20] : 
{d^2_y2 , dxy, d^2, + Pxl)lV^, Vv^} and 

{dxz, dyz, {p^x\ - Px}y)l'j2, {p^z'’ + P^z'’)/ V2}- 

The full group of the point symmetries of M0S2 SLVD 
with our considered VDs is D^h = Csy (8) Z2. Thus the 
states bound to the VDs, even and odd with respect to 
cT/i, must transform according to Ai^2 and the irre¬ 
ducible representations of C^y. Respectively, the bound 
states must appear as singlets and doublets. It should be 
noted that such classification holds if the overlap between 
states bound to different VDs is absent. 



FIG. 2. Bandstructure of different kinds of VDs. The Fermi 
level is set at cf = 0 eV. Red (blue) lines show odd (even) 
states. Arrows indicate transitions corresponding to reso¬ 
nances in XF(ct;) (see Eq. @ and below) shown in Fig. 

(a) Mo-vacancy, (b) S 2 -vacancy. Here singlets and doublets 
are in different half-planes with respect to the Fermi level. 
Therefore there are no tt transitions, (c) S 2 -vacancy with 
charge added to M 0 S 2 SLVD, raising the Fermi level such 
that TT transitions become allowed, (d) 3xMoS2 QAD. 


The simplest model describing the arrangement of the 
bound states is the TBM considering only the atoms on 
the edge of the VD, which is inferred by the small local¬ 
ization radius of the bound states. For the case of the 
hexagonal 3xMoS2 QAD one finds 

e = e ± y/Se‘^ + 4 |tp cos^(^), ( 1 ) 


where e = (cmo + ^s) /2, ( 5 e = (cmo - es) /2 and ^ = 0 
for singlet states (invariant with respect to C3 rotations) 
and ^ = 27 r /3 for doublets (states aquiring the phase 
factor exp(±i 27 r/ 3 )). Here cmo, are the phenomeno¬ 
logical parameters describing the energy of the electron 
on Mo and S atoms, respectively, and t is the hopping 
parameter. Equation 0 correctly reproduces the se¬ 
quence {A[^E\E',A[} for even and {A'{,E"^E"^A'{} 
for odd states, i.e. {singlet, doublet, doublet, singlet}, 
while traversing the gap from the bottom of the con¬ 
duction band down over the bound states (see Fig. It). 

It may appear that this model contradicts the numer¬ 
ical results for the Mo-vacancy, where the numerical cal¬ 
culations show only 5 bound states (see Fig. [^). How¬ 
ever, the TBM model suggests that in addition to the 
bound states appearing inside the gap Eg\\ of M0S2 SL 
there must be states, in this particular case a singlet 
state, hidden inside the bands. Indeed, there is such 
a state inside the valence band at energy e ^ — 0.5 eV 
below the top of the valence band (see Fig. IT)- 

Moreover, the parameters cmo and es are deter¬ 
mined by the microscopic Hamiltonian, e.g. es = 




n 


k(S) 


y. Thus, we can expect that there is a va¬ 
riety of states bound to VDs besides the ones inside the 
gaps E^||, Eg_i. An example of such states is provided 
by the case when the bound state is made of Sulfur’s 
5-orbitals (see Fig. [^) at energy ^ —12 eV below the 
Fermi level. 

Dirac model The TBM considered above relates the 
structure of the spectrum of the bound states to the sym¬ 
metry of the VD. Due to the fact that its parameters 
should be fitted to the energies of electron states obtained 
by other means, however, it cannot explain neither the 
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FIG. 3. Examples of the electron probability distributions in 
states bound to Mo-vacancy. Single super-cell is shown, black 
and gray dots indicate positions of Mo and S atoms, respec¬ 
tively. The distances are measured in A. (a) Odd doublet at 
energy 0.5 eV above the top of the valence band, (b) odd 
singlet inside the valence band at energy 0.5 eV below the 
top of the valence band, (c) deep defect state at energy ^ —12 
eV below the Fermi level formed by Sulfur’s s orbitals. 


smallness of the localization radius of the bound states 
nor their energies. In particular, it cannot explain why 
the odd states may form bound states inside the gap . 
These features, however, can be understood with the help 
of an analysis of circularly symmetric QADs based on 
the Dirac equation which emerges as the two-band model 
within the k • p-approximation near the iF-point of the 
Brillouin zone of M0S2 SL. 

Considering two bands with the energy separation 2A 
between them, the equation describing the spatial dis¬ 
tribution of the pseudo-spin has the form Hr^r = 
where r = ±1 enumerates the valleys, the energy refer¬ 
ence level is chosen to be positioned at the center be¬ 
tween the bands, and 1 -Lt = tcJz^ + va • p. Assum¬ 
ing that the QAD has circular shape, we rewrite this 
equation in the polar coordinates and withjrespect to 
= exp(icr;2 0/2)$T-Y^ we obtain where 


Hr = TazA — iv 


d Id 
dr ^ dcj) 


( 2 ) 


The solution is subject to the condition of vanishing ra¬ 
dial component of the probability current at the bound¬ 
ary [ 2 T 1 [ 22 ] (^(ro) \ nB ' cr\ ^(ro)) = 0, where tq is the ra¬ 
dius of the QAD and is the unit vector perpendicular 
to the boundary. The straightforward implementation of 
such boundary condition is provided by the infinite mass 
model [ 23 ], where the QAD is represented by a region 
with renormalized width of the gap A ^ A(1 + d{r))^ 
with d{r) = 0 for r > ro and d{r) oc when r < ro- 
In our case we identify A = Within this model 

the boundary condition is satisfied if |^(ro)) oc \ry), 
i.e. the pseudo-spin |^(ro)) is tangent to the bound¬ 


ary of the QAD. Next, observing that az \ry) = \—ry) 
and ax \ry) = ir \ —ry), one can see that Hr has an an¬ 
gularly independent solution |^(ro)) oc exp(—r/rc) \ry) 
corresponding to e = 0, which exponentially decays for 
r > ro with the localization radius Vc = A/v. 

Thus, independently of its radius the QAD may sup¬ 
port a bound state with very short localization length and 
with the energy in the middle between the energies of the 
coupled bands. Comparing this finding to the distribu¬ 
tion of energies of even and odd bands [20] the conclusion 
can be drawn that, indeed, both even and odd bands may 
support bound states with the energy near the energy of 
the Fermi level, i.e. inside the gap E^y of M0S2 SL. 

Optical spectrum. In view of nontriviality of the ap¬ 
pearance of odd bond states inside the gap E^y it is im¬ 
portant to note that the presence of the bound states of 
different parities manifests itself in the optical spectrum 
of M0S2 SLVD. Therefore, they are available for a direct 
experimental observation. 

When VDs form a superlattice the problem of the op¬ 
tical response can be approached along the same line as 
for single layered systems [ 24 ]. Let k be a point in the 
first Brillouin zone of the superlattice. At this point the 
electron wave function satisfies 

e„(k) IV’n(k)) = ■H(k) |V’n(k)), ( 3 ) 

where n enumerates the super lattice bands. Implement¬ 
ing the k • ^-approximation of Eq. in the usual way 
and using the Peierls substitution we obtain the Hamilto¬ 
nian of interaction with the elecromagnetic field Hem = 
• p. Treating Hem as a perturbation within the lin¬ 
ear response theory we find the Kubo-Greenwood optical 
susceptibility (see e.g. Ref. [ 25 ]) 

2kA r 

- 2 E / c^kP„,„/(k) (8)P„/,„(k) 

n,n' 

^ /[en(k)] -/[en/(k)] 

e„(k) - e„/(k) - u - 

where /(e) is the Fermi distribution, (g) denotes the tensor 
product and Pn,n'(k) = {2pn{^)\p |'0n'(k)). 

The appearance of the states inside the gap E^y of 
M0S2 SL leads to resonances at frequencies of corre¬ 
sponding transitions. Several transitions, however, are 
prohibited due to symmetry, i.e. when Pn,n' (k) does not 
transform according to the symmetric representation of 
the symmetry group of the superlattice. Pn,n' (k) trans¬ 
forms according to /(Cs^)^ (g) ^(^2)^ (g) /(P), where I{G) 
and /(P) denote irreducible representations of group G 
and the momentum operator P, respectively, and pow¬ 
ers are shorthand notations for direct products. One 
needs to consider separately the in-plane and out-of-plane 
components of Pn,n' because they transform according 
to different irreducible representations of namely, 

E' = E 0 A' and A2 = Ai (g) A", respectively. Taking 
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FIG. 4. (a) Resonances of Inixy (cc;) (blue) and Inix±(ct;) (red) 
in pristine M 0 S 2 SL. (b) M 0 S 2 SLVD with Mo-vacancy: T 3 
resonance in x±{^) is due to the odd singlet state inside the 
valence band, (c) M 0 S 2 SLVD with S 2 -vacancy: If uncharged 
(green line), tt transitions are suppressed; if charged (red line), 
TT transitions are allowed, (d) M 0 S 2 SLVD with 3 xMoS 2 . 


into account the multiplication rules for Csy'. 2 — 

0 E = E and = Mi 0 M 2 0 E, we find that the 
out-of-plane component of Pn,n'(k), which gives rise to tt 
transitions, is nonzero only between odd and even states 
of the same multiplicity (either between singlets or be¬ 
tween doublets), while the in-plane components, which 
lead to a transitions, are nonzero for all states of the 
same parity. Thus x is diagonal in the basis spanned by 
ex,y,z and is isotropic in the plane of the layer and, thus, 
is characterized fully by two eigenvalues X|| and x±- 

The numerical results for the optical spectrum are 
shown in Fig. The difference between X-l(<^) pris¬ 
tine M 0 S 2 SL and for M 0 S 2 SLVD is drastic. For pristine 
M 0 S 2 SL the lowest energy transition Tq yielding nonzero 
Xx (see Fig. Ef) corresponds to the transition between 
the top of the valence band to the CB 0 1 band with 
energy 3.2 eV (see Fig. [^). In turn, for M 0 S 2 SLVD the 
lowest energy resonance is due to the transition between 
bound states of the same degeneracy with the energy dif¬ 
ference smaller than 1 eV (see Fig. g), c, and d). This 
result is in stark contrast to true 2D systems where tt 
transitions are absent, and the effect of non-zero small 
thickness may be expected to be observed at energies at 
least significantly higher than those characteristic to a 
transitions. In addition, the selection rules governing tt 
transitions present a great opportunity for experimental 
characterization of states bound to VDs. 

The qualitative picture based on the symmetry prop¬ 
erties establishes the connection between the main fea¬ 


tures of spectrum of the bound states and the optical 
response. For example, for the Mo-vacancy the T 3 res¬ 
onance in X-l(<^) involves a bound state hidden in the 
valence band. In the case of S 2 -vacancy the symmetry 
analysis predicts that x_l(co’) is featureless at low energies 
due to the smallness of // — /i, which is confirmed by the 
numerical calculations (see Fig. g)- The Fermi level 
however, can be shifted by means of a gate voltage to 
lie between equally degenerate states with different par¬ 
ities, which support transitions contributing to x±{^)- 
Then x± (<^) should demonstrate a low-energy resonance. 
In the numerical simulations we modified the position of 
the Fermi level by adding charge to the whole layer by 
means of a charge concentration of 1.4 x 10^^cm“^, lead¬ 
ing to a resonance in x±{^) (see Fig. g)- Fig. 1^ shows 
the resonances due to transitions between states bound 
to 3 XM 0 S 2 . 

Considering the case when there is no overlap be¬ 
tween bound states of neighboring VDs, we can use 
ImX±,Ti = p{d‘^Jeoh)^/[y‘^^{en-en' for the transi¬ 

tion Ti of a dilute gas of VDs [26], where p = N/V and N 
are the concentration and number of VDs, respectively. 


dTi = (r) I r (i*)^ deuotes the dipole moment of 


the transition T^. This formula is in excellent agreement 
with the numerical calculations shown in Fig. [^ for su¬ 
percell sizes from 7x7 up to 13x13. The peak for the 5x5 
supercell does not follow this formula because the over¬ 
lap between neighboring VDs is substantial, which leads 
to a peak shift and homogeneous peak broadening due to 
the formation of minibands. Additional inhomogeneous 
peak broadening is expected due to random distribution 
of VDs in the M 0 S 2 SL. 



FIG. 5. (a) Magnitude of x^(^) ^ function of VD concen¬ 

tration p for Mo-vacancy around the peak T 2 (see Fig. 1^) 
for the five supercells 5x5, 7x7, 9x9, 11x11, and 13x13. (b) 
Linear dependence of x^(^) ^it maximum of peak T 2 . 

Conclusion. We show that in order to describe the 
electron states bound to VDs in M 0 S 2 SLVD, it is neces¬ 
sary to consider odd states, which lead to the appearance 
of resonances in the out-of-plane optical response x±(<^)- 
Our results pave the way to the optical characterization 
of VDs in TMDC SLVD, which is of utmost importance 
for the future realization of high-performance electronic 
and optoelectronic devices based on TMDC SLs. 
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